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Spin Polarization at Semiconductor Point Contacts in 
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Semiconductor point contacts can be a useful tool for producing spin-polarized currents in 
the presence of spin-orbit (SO) interaction. Neither magnetic fields nor magnetic materials 
are required. By numerical studies, we show that (i) the conductance is quantized in units of 
2e 2 /h unless the SO interaction is too strong, (ii) the current is spin-polarized in the transverse 
direction, and (iii) a spin polarization of more than 50% can be realized with experimentally 
accessible values of the SO interaction strength. The spin-polarization ratio is determined by 
the adiabaticity of the transition between subbands of different spins during the transport 
through the point contacts. 

KEYWORDS: point contact, Rashba, spin-orbit interaction, conductance quantization, spin filter, Landau- 
Zener 



Generating spin-polarized currents in semiconductors 
is an important issue for the development of spin-based 
electronics, "spintronics." 1 To manipulate electron spins, 
the Rashba spin-orbit (SO) interaction is useful since 
its strength is locally controllable by applying an elec- 
tric field. 2 ' 3 Several spin- filtering devices for producing 
the spin currents have been proposed utilizing the SO 
interaction, e.g., three-terminal devices related to the 
spin Hall effect, 4,5 a triple-barrier tunnel diode, 6 a one- 
dimensional system with a magnetic field, 7 and a three- 
terminal device for the Stern-Gerlach experiment using 
a nonuniform SO interaction. 8 

In the present paper, we theoretically study the bal- 
listic transport through a semiconductor point contact 
(quantum wire with a narrow constriction) in the pres- 
ence of Rashba SO interaction. In its absence, it is 
well known that the conductance is quantized in units 
of 2e 2 /h when the constriction changes gradually in 
space. 9,10 By numerical studies, we show that the con- 
ductance is quantized even with the SO interaction and 
that the current is spin-polarized in the transverse direc- 
tion, in the absence of magnetic field. We demonstrate 
that the polarization ratio can be more than 50% in 
InGaAs heterostructures. The spin current is obtained 
generally, e.g., in GaAs heterostructures, if the condi- 
tion later described by eq. (9) is fulfilled with not only 
Rashba but also Dresselhaus SO interaction. 11 As spin 
filters, two-terminal devices with point contacts are easy 
to fabricate on semiconductors, compared with other de- 
vices that have been proposed. 4-8 

We consider a two-dimensional electron gas confined 
in the z direction. The electric field in the z direction 
results in the Rashba SO interaction, 



H 



RSO 



■{Py&x -Px&y), 



(1) 



where a x and a y are Pauli matrices. (The Dresselhaus 
SO interaction is discussed later.) We use a dimensionless 
parameter, k a /kp, where 



k a = ma/h 2 



(2) 



with to being the effective mass and hp is the Fermi 
wavenumber. 12 In InGaAs heterostructures, a — (3 ~ 
4) x KT 11 eVm and A R = 2k F a = 15-20 meV 
[k a /kp = Ar/(4E f ) w 0.1]. 13,14 Electron-electron in- 
teraction and impurity scattering are neglected. Elec- 
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trons propagate in a quantum wire along the x direc- 
tion, with width Wo in the y direction. A hard-wall 
confinement potential U(y) is assumed, U(y) = for 
—Wo/2 < y < Wo/2 and oo otherwise. For a narrow 
constriction around x = y = 0, we consider an extra 
potential at —L\ < x < L2, which is analogous to that 
adopted in ref. 15: 



V(x,y) 



1 



A 



L, 



(±[y-y ± (x)}), (3) 



with 



y±(x) = ± 



W 



1 — cos ■ 



L x 



(4) 



where L x — L x (— £1 < x < 0) and L x — L 2 (0 < x < 
L 2 ). 6{t) is a step function \6{t) — 1 for t > 0, for 
t < 0]. For fixed x, V(x,y) is flat at y~(x) < y < y+(x) 
and grows with increasing \y — y±(x)\ at y > y+{x) or 
y < y~(x). On a line of y = 0, the potential height is 
given by (Vb/2)[l+cos(7ra;/L 2 ;)], being maximal at x = 0. 
In the present paper, we fix Wo = 4Af and A = Ap, 
where Af is the Fermi wavelength (Af = 2n/kF). 

Numerical calculations are performed using a tight- 
binding model on a square lattice (— L\ < x < L2, 
-Wo/2 < y < Wo/2), following refs. 15 and 16. 17 The 
transmission coefficients are evaluated for incident elec- 
trons from the left side of the constriction (x < — L\) 
to the right (L2 < x), using the Green function's recur- 
sion method. 15 They yield the conductance G through 
the Landauer formula. 

Figure 1 shows the calculated results when L\ = L2 = 
4Af- In Fig. 1(a), we plot G as a function of Vq, poten- 
tial height at x = y — 0. The conductance quantization 
is clearly observed when k a /k-p — 0.25 (solid line) as 
well as in the absence of SO interaction (broken line). 
The conductance quantization is broken for k a /k-p > 0.5 
(not shown here), which is addressed later. We divide 
the output current into two components, one carried by 
spin-up electrons in the y direction (S y = 1/2) and the 
other by spin-down electrons (S y = —1/2). Figure 1(b) 
presents each conductance, G±, as a function of k a /k?. 
Vq = 0.7-Ef where the total conductance is at the first 
plateau (G = G+ + G_ = 2e 2 /h). The spin polarization 
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Fig. 1. Numerical results of the conductance G through a point 
contact with L\ = L2 = 4Ap. (a) G as a function of Vq, po- 
tential height at x = y = 0. The strength of the spin-orbit 
(SO) interaction is k a /kp = 0.25 (solid line). The broken line 
indicates G in the absence of SO interaction. Inset: Schematic 
drawing of our model (Wo = 4Ap). (b) The conductance G± 
for electrons with S y = ±1/2 in the output current, as a func- 
tion of k a /kp. Vo = 0.7 Ep where G is at the first plateau 
(G = G+ + G_ =2e 2 /h). 



in the y direction, (G+ — G_)/(G-|_ + G_), increases with 
an increase in k a /kp. Note that incident electrons are un- 
polarized in eight conduction channels per spin direction. 
These results indicate that (i) the point contact works as 
a spin filter, (ii) the polarization ratio is about 30% with 
experimental values of k a /kp « 0.1 in this case, and (iii) 
the conductance is still quantized. 

These calculated results can be understood as follows. 
We divide the Hamiltonian into two parts: H = Hq + H', 



Ho = 



H' 



1 

2m 

Q 



(Pl+P 2 y) ~ 



-PyOx- 



-p x <J y + V(x,y) + U(y),(5) 
(6) 



We treat H' as a perturbation. When V(x, y) is indepen- 
dent of x, the eigenstates and eigenvalues of Hq are given 
by 



1pn,k,± 
E n ,±{k) 



Akx 



i Pn(y)x±, 



n 2 

2m 



(7) 
(8) 



respectively, where <p n (y) are eigenstates of Py/(2m) + 
V(y) + U{y) with eigenvalues e n (n = 1,2,3, •••). x± 
are eigenstates of a y , representing spin-up or -down 
states. The dispersion relations of E n ^±(k) (subbands) 
are schematically shown in Figs. 2(a) and 2(b). The 
Fermi level is denoted by horizontal line A in the quan- 
tum wire outside of the constriction, where V(x,y) = 0. 
(The number of channels is three per spin direction in 
the figures, which is smaller than that in our numerical 




Fig. 2. Schematic drawings of subbands, E n ^(k) (solid lines) 
and E n , — (k) (broken lines) with n = 1,2,3, in a quantum wire 
with a component of SO interaction, —(a/H)p x a y . All the cross- 
ings between E ny —(k) and E n i +(fe) (n < n') take place above 
the Fermi energy E-p in case (a), whereas some of them appear 
below Bp in case (b). Horizontal lines indicate E-p (relative to 
the subbands) corresponding to three positions of electrons in 
the point contact, shown in the inset in (a), (c) A vicinity of the 
crossing between i?i _(fc) and i?2,+(&), which is surrounded by 
a square in (b). The subbands are mixed with each other by the 
other component of SO interaction, (a/h)p y cr x . 



studies.) There are two situations regarding the cross- 
ings between E n ^(k) and E n >. + (k) (n < n'). When 
k a /kY < (3/16)(Af/Wo) 2 , all the crossings take place 
above E F [Fig. 2(a)]. When 



k a /k F > (3/16)(A F /Wo) 2 



(9) 



some of them appear below Ep [Fig. 2(b)]. 18 

Now we consider the transport of electrons through 
a constriction when the conductance shows the first 
plateau. Except in the vicinities of the above-mentioned 
crossings, we assume an adiabatic transport in which 
the wavefunction ^p n ,k,± changes gradually remaining the 
quantum numbers of transverse motion n and of spin 
zb. ' The wavenumber k changes with x, which is de- 
termined by an intersection between the subband and 
Ep; positive k's for incident electrons. Let us begin with 
the case of Fig. 2(a). As electrons propagate from the 
wire to a narrow constriction (x < 0), the subbands shift 
upwards. Alternatively, we move E-p downwards in Fig. 
2(a). Similarly, at x > 0, we move E-p upwards. [Sep- 
arations between the subbands increase (decrease) with 
x at x < (x > 0), which is not shown in the figure.] 
Before the injection into the constriction (position A), 
there are six conduction modes, (1, ±), (2, ±) and (3, ±). 
At position B, modes (1,±) and (2,±) are conducting, 
whereas modes (3, ±) have been completely reflected. At 
the narrowest region (position C), only modes (1, ±) ex- 
ist. At x > 0, the modes propagate with the transmission 
probability of unity, which results in the conductance 
quantization, G = 2e 2 /h. The small perturbation of H' 
does not play an important role. No spin polarization is 
observed in this case. 

In the case of Fig. 2(b), the situation is different. Elec- 
trons pass by the crossings twice, once at x < and once 
at x > 0. Let us look at the crossing between modes 
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(a) 



Fig. 3. Conductance G± for electrons with S y = ±1/2 in the 
output current, as a function of k a /kp (solid lines for G+ and 
broken lines for G_). Vq = 0.7 E-p in eq. (3) where G is at the 
first plateau (G = G+ + G_ = 2e 2 /h). L\ = L 2 = 4A F . A term 
of (a/ h)p y cr x in -Hrso i s taken into account (a) only at x > 
or (b) only at x < 0, whereas the other term, —(a/h)p x cr y , is 
considered in the whole region. The thin lines indicate G± when 
both the terms are present in the whole region. 



(1, — ) and (2, +). These subbands are mixed by the per- 
turbation H', as shown in Fig. 2(c). Hence these modes 
change to each other with a transition probability P 
when electrons pass through the crossing. Around the 
first pass (x < 0), both modes are occupied by electrons 
just before the modes cross. Then spin-up electrons are 
flipped to spin-down with probability P, while spin-down 
electrons are flipped to spin-up with the same probabil- 
ity. Accordingly, no spin polarization takes place. Around 
the second pass (x > 0), on the other hand, mode (2, +) 
is empty while mode (1, — ) is full of electrons just before 
the modes cross. Then spin-down electrons in the lat- 
ter mode are spin-flipped to spin-up in the former mode 
with probability P. The spin-up electrons in mode (1, +) 
are transmitted through the constriction without pass- 
ing by any mode crossing. Consequently we obtain the 
spin-polarization ratio of [(1 + P) — (1 — P)]/2 = P. 

The transition probability P is evaluated using the 
Landau-Zener theory: 19,20 



P = 1 - exp(-27rA), 



(10) 



where A = J 2 /(H\v\) represents the degree of adiabatic- 
ity J = |(2, +|ir / |l, — )| and?, = £[E(2,+) - E(l, -)}, 
the velocity of the change of level spacing. P = 1 for 
A = oo in the adiabatic limit, whereas P — for A = 
in the sudden-change limit. If V(x, y) is a hard-wall po- 
tential with width W(x) in the y direction, instead of by 
eq. (3), A is estimated as 



A 



1 dW 
W~dx 



(11) 



apart from a numerical factor. W and its derivative 
should be evaluated at x where electrons pass by 
the crossing. In more gradual point contacts (smaller 
e l ec t rons P ass through the crossing more adi- 
abatically (larger A). Then the larger spin- flip probabil- 
ity P, and thus the larger spin-polarization ratio, is ex- 
pected. 

In spite of the spin polarization, the total conductance 
is not affected by the mode crossings, unless the SO in- 
teraction is too strong. In Fig. 2(b), both the derivatives 
of _Ei _(fc) and E 2 ,+(k) are positive at their crossing. 
Since the group velocities have the same sign, the tran- 
sition from one mode to the other is not accompanied 
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Fig. 4. Conductance G± for electrons with S y = ±1/2 in the 
output current, as a function of k a /kp (solid lines for G+ and 
broken lines for G— ). Vb = 0.7-Bf in eq. (3) where G is at the 
first plateau (G = G+ + G_ = 2e 2 /h). (a) L\ is fixed at 4A F , 
whereas Li is changed as (A) 4Ap, (B) 8Ap and (C) 12Ap. (b) 
Li is fixed at 4Ap, whereas L\ is changed as (A) 4Ap, (B) 8Ap 
and (C) 12A F . 



by a reflection (forward scattering). If the SO interac- 
tion were so strong that the crossing occurred at k < k a , 
backward scattering could take place, which would de- 
stroy the conductance quantization. The condition that 
the backward scattering does not take place is given by 



2h I k I l /m < e 2 - ei, or 



k a /k F < (V3/4)(A F /W0 



(12) 



for the hard-wall confinement of width W(x). In eq. (12), 
W is the width of confinement where modes (1,—) and 
(2,+) cross (Af <W< Wo). [In our numerical studies, 
this condition seems to be satisfied with k a /kp < 0.25. 
We could still observe a spin-polarized current when eq. 
(12) does not hold, as discussed later.] 

We have so far presented a simple theory to explain 
the numerical results in Fig. 1, on the basis of the adia- 
batic approximation for the transport through the point 
contact and perturbative treatment of H' in eq. (6). To 
verify this theory, we further perform numerical studies 
and present the results in Figs. 3 and 4. 

In our theory, the spin-flip by H' is important for the 
spin polarization during the transport from a narrow re- 
gion to a wide region [x > 0). To confirm this idea, we 
examine two situations in Fig. 3: (a) H' is present only 
at x > or (b) only at x < 0. Indeed we observe the spin 
polarization in situation (a), whereas the polarization is 
not visible in situation (b). 

The adiabaticity of the spin-flip transition can be con- 
trolled by changing the shape of point contacts, e.g., 
(Li,L 2 ) in our model. With increasing L 2 , \^^-\ de- 
creases in eq. (11) and hence P in eq. (10) increases. 
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As a result, a larger spin-polarization is expected. Fig- 
ure 4(a) shows the numerical results with (A) (L\, L2) = 
(4A F ,4A F ), (B) (4A F ,8A F ), and (C) (4A F ,12A F ). The po- 
larization ratio increases with an increase in L 2 , in accor- 
dance with our theory. In case (C), we observe a polar- 
ization of 60% with experimental values of a in InGaAs 
hctcrostructurcs. In Fig. 4(b). we change Li with a fixed 
Li. The spin polarization is not influenced by L\, which 
is consistent with the previous discussion. 

In conclusion, we have examined the ballistic trans- 
port through semiconductor point contacts in the pres- 
ence of Rashba SO interaction. We have observed a spin- 
polarized current although the conductance is still quan- 
tized. The spin-polarization ratio is determined by the 
adiabaticity of the transition between subbands of dif- 
ferent spins, which is characterized by eq. (11), during 
the transport from a narrow region to a wide region. 

We have examined a quantum wire of width Wq = 4A F 
with a narrow constriction and demonstrated that the 
polarization ratio can be 60% in InGaAs heterostruc- 
tures. Generally, a condition to generate the spin current 
is given by eq. (9). In GaAs heterostructures, Wo > 6A F 
is required with a = 0.05 x 10~ n cVm (k a /kp ~ 
0.005). 21 Our mechanism also works with Dresselhaus 
SO interaction 



-ffDSO = "^T" ( Vx®x + Py cr y)- 



(13) 



Then the output current is spin-polarized in the x di- 
rection. If i?R,so and i?Dso coexist, 21 the eigenstates of 
cujy + a'cr x determine the direction of spin polarization. 

We make a few remarks, (i) The spin filtering effect is 
expected at the higher plateaus of conductance as well as 
at the first plateau. At the second plateau, for example, 
the crossings between modes (1, — ) and (n, +) and those 
between (2, — ) and (n, +) (n > 2) work for the spin po- 
larization if they appear below Ep. Note that the cross- 
ing between modes (1, — ) and (2, +) docs not work since 
both modes are occupied just before the modes cross even 
at the second time, (ii) A steplike structure of the spin 
polarization is seen as a function of a in Figs. 1(b) and 4. 
This reflects the number of crossings between (1, — ) and 
(n, +) with n > 1 below Ep. With an increase in a, more 
crossings appear below Ep, which increases the spin- flip 
probability, (iii) In our model, the conductance seems to 
be influenced by the backward scattering for k a /kp > 0.5 
[eq. (12) does not hold]. Around k a /kp = 0.7, we observe 
the conductance plateaus with a small fluctuation, which 
might be due to some resonant states around the point 
contact with backward scattering. A spin current is still 
obtained with a polarization ratio of ~90%. When the 
SO interaction is increased further, the plateaus begin to 
break. The investigation of this regime in more detail is 
beyond the scope of this paper. 

Finally, we discuss the observation of the spin- 
polarized current produced by our mechanism. The spin 
polarization has to be directly measured since the con- 



ductance is not influenced. Possible experiments are an 
injection of the spin current into dilute magnetic semi- 
conductors, an injection into a spin detector, 22 ' 23 or an 
optical measurement. 24 An indirect measurement may 
be available owing to the fact that mode (1,—) is con- 
verted to upper modes, (n, +) with n > 1, in the spin- flip 
processes. If the obtained current is injected into another 
point contact in the absence of SO interaction, the height 
of the first plateau is suppressed to ~ (2 — P)e 2 /h when 
the polarization ratio is P. 

The authors gratefully acknowledge discussions with 
D. Matsubayashi, H. Yokouchi, J. Yamauchi, Y. Tokura 
and G. E. W. Bauer. This work was partially supported 
by a Grant-in- Aid for Scientific Research in Priority Ar- 
eas "Semiconductor Nanospintronics" (No. 14076216) of 
the Ministry of Education, Culture, Sports, Science and 
Technology, Japan. 



13 
14 

15 
16 
17 

18 



19 
20 
21 



22 
23 

24 



I. Zutic, J. Fabian and S. Das Sarma: Rev. Mod. Phys. 76 
(2004) 323. 

E. I. Rashba: Fiz. Tverd. Tela (Lenningrad) 2 (1960) 1224 
[Solid State Ion. 2 (I960) 1109]. 

J. Nitta, T. Akazaki, H. Takayanagi and T. Enoki: Phys. Rev. 
Lett. 78 (1997) 1335. 

A. A. Kisclev and K. W. Kim: Appl. Phys. Lett. 78 (2001) 
775. 

T. P. Pareek: Phys. Rev. Lett. 92 (2004) 76601. 

T. Koga, J. Nitta, H. Takayanagi and S. Datta: Phys. Rev. 

Lett. 88 (2002) 126601. 

Y. Stfcda and P. Seba: Phys. Rev. Lett. 90 (2003) 256601. 
J. Ohe, M. Yamamoto, T. Ohtsuki and J. Nitta: cond- 
mat/0409161. 

B. J. van Wees, H. van Houten, C. W. J. Beenakker, J. G. 
Williamson, L. P. Kouwenhoven, D. van der Marel and C. T. 
Foxon: Phys. Rev. Lett. 60 (1988) 848. 

A. Kawabata: J. Phys. Soc. Jpn. 58 (1989) 372. 

G. Dresselhaus: Phys. Rev. 100 (1955) 580. 

fcp (Ep) is the Fermi wavenumber (energy) in the absence 

of SO interaction; Ep = h 2 kp/(2m). In its presence, we shift 

the Fermi energy by —h 2 k 2 t /(2m) to keep the electron density 

constant. 

D. Grundler: Phys. Rev. Lett. 84 (2000) 6074. 

Y. Sato, T. Kita, S. Gozu and S. Yamada: J. Appl. Phys. 89 

(2001) 8017. 

T. Ando: Phys. Rev. B 44 (1991) 8017. 
T. Ando and H. Tamura: Phys. Rev. B 46 (1992) 2332. 
The lattice constant a is chosen as Xp/a = 8 when k a /kp < 
0.25. A smaller value is required for a when k a /kp > 0.25. 
Equation (9) is an approximate expression that is valid for 



fcjl > {2n/W f 



It is 



■J 



exactly given by k a [^J k F — (tt/Wq) 2 



kp - 



(27r/W )2]/2>(3/4)(7r/^ ) 2 . 
L. Landau: Phys. Z. Sowjetunion 2 (1932) 46. 
C. Zcner: Proc. R. Soc. London, Ser. A 137 (1932) 696. 
J. B. Miller, D. M. Zumbiihl, C. M. Marcus, Y. B. Lyanda- 
Geller, D. Goldhaber-Gordon, K. Campman and A. C. Gos- 
sard: Phys. Rev. Lett. 90 (2003) 76807. 

J. A. Folk, R. M. Potok, C. M. Marcus and V. Umansky: 
Science 299 (2003) 679. 

J. M. Elzerman, R. Hanson, L. H. Willcms van Beveren, B. 
Witkamp, L. M. K. Vandcrsypen and L. P. Kouwenhoven: 
Nature 430 (2004) 431. 

Y. K. Kato, R. C. Myers, A. C. Gossard and D. D. Awschalom: 
Nature 427 (2003) 50. 



